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The topological charge is studied on lattices of large physical volume and fine lattice spacing. We 
illustrate how a parity transformation on the SU(3) link- variables of lattice gauge configurations 
reverses the sign of the topological charge and leaves the action invariant. Random applications of 
the parity transformation are proposed to traverse from one topological charge sign to the other. The 
transformation provides an improved unbiased estimator of the ensemble average and is essential in 
improving the ergodicity of the Markov chain process. 
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I. INTRODUCTION 

The Markov chain process is the corner stone of mod- 
ern importance-sampling techniques for generating field 
configurations in the numerical simulation of quantum 
field theories. Central to the process is ergodicity; the 
ability to move through configuration space from any 
particular field configuration to any other field configura- 
tion with finite probability. Vacuum expectation values 
of observable operators are estimated from the average of 
field configurations selected with a probability given by 
the exponentiation of the Euclidean action, exp(— Se), 
governing the quantum field theory. 

It is essential that these field configurations are selected 
based on the action alone, unbiased by the field configu- 
ration used to initiate the Markov chain. In practice this 
is done by updating the field configurations for thousands 
of sweeps through the lattice, monitoring the autocorrela- 
tion of various observables, and selecting a new represen- 
tative field configuration only after significant evolution 
through configuration space. 

The topological charge of a gauge configuration in lat- 
tice QCD has already been identified as a particularly dif- 
ficult quantity to evolve in the Markov chain process, dis- 
playing unusually long autocorrelation times 0, 0, S Q • 
Some gauge actions such as the renormalization-group 
block-transform based DBW2 action Q are notorious 
for locking in the topological charge at fine lattice spac- 
ings H, 0| . Difficulties associated with the Iwasaki gauge 
action are presented here. 

In this study we address an aspect of the topological 



charge evolution problem that will present itself for any 
gauge field action. In particular, we address a problem 
associated with the approach to the infinite-volume con- 
tinuum limit in numerical supercomputer simulations of 
SU(3) gauge theory and QCD in general. 

As the physical volume, V, of the lattice increases, the 
average value of the squared topological charge increases, 
in accord with the topological susceptibility, given by 
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In quenched QCD, \ is related to physical hadron masses 
via the large N c Witten-Veneziano relation jjj 
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In full QCD, the quark mass dependence of the topolog- 
ical susceptibility may be related to pseudoscalar meson 
propert ies via the Gell-Mann-Oakes-Renner relation as 
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On a sufficiently large volume lattice, the distribution of 
the topological charge is expected to be Gaussian. Re- 
gions of non-trivial topological charge density are uncor- 
related for sufficiently large separations and a normal dis- 
tribution will result. The distribution of Q for the Wil- 
son gauge action has been found to be Gaussian to a very 
good approximation for both quenched QCD and 
full QCD where the Wilson gauge action is complemented 
by the Wilson fermion action [l( or the Wilson-clover ac- 
. , . tion 0. 
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FIG. 1: Histogram illustrating the number of configurations in 
our ensemble of 250 configurations having a particular topo- 
logical charge, Q. A suppression of Q ~ configurations is 
observed. 
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FIG. 2: The time evolution of the topological charge plotted 
as a function of simulation time as represented by the config- 
uration number. 1000 pseudo-heatbath sweeps are performed 
between each gauge field sample. 



Figure ^ displays a histogram of the topological charge 
from an ensemble of 250 quenched Iwasaki [7J, |8( gauge 
configurations at (3 = 9.1674 on a 28 3 x 44 lattice, where 
the lattice spacing a is 0.113 fm. The topological charge is 
calculated using the highly-improved, three-loop 0{a A )- 
improved lattice field-strength tensor on 10-sweep 
cooled configurations obtained with a three-loop 0(a A )- 
improved action. The gauge configurations represented 
in Fig. ^are separated by 1000 pseudo-heatbath sweeps in 
an attempt to obtain uncorrelated configurations. Fig.|2 
illustrates the time evolution of the topological charge 
plotted as a function of simulation time, represented by 
the configuration number. Acceptable movement of the 
topological charge is indicated by the rapidity and am- 
plitude of the oscillations. To the best of our knowl- 
edge, this is the first time a double-peaked structure 
in the probability distribution of the topological charge 
has been revealed for a renormalization-group improved 
gauge action. It would be interesting to examine this 
distribution for the DBW2 action. 

For any lattice action, the distribution of the topo- 
logical charge will broaden as the volume of the lattice 
increases, such that 

(Q 2 )=xv. (4) 

As the infinite- volume continuum limit is approached and 
correlation lengths diverge, it will become increasingly 
difficult to evolve the topological charge over the broad 
range demanded by Eq. In particular, a symmet- 

ric distribution about Q = is required to preserve the 
symmetries of QCD. 

Recognizing the need for improvement in the evolution 
of the topological charge in the Markov chain process, 
we present a simple transformation in the following sec- 
tion, that may be applied to a gauge-field configuration of 
quenched QCD or dynamical-fermion QCD that changes 



the sign of Q while leaving the action invariant. The 
transformation provides an improved unbiased estimator 
of the ensemble average and is essential in improving the 
ergodicity of the Markov chain process. 

II. GAUGE-LINK TRANSFORMATION 

The index theorem relates the topological charge 
Q to the chirality index of the Dirac operator D on a 
continuum 4-torus as 

Q = index(_D) = n_ — n+ . (5) 

Here n+ and n_ are the number of exact zero eigen- 
modes, D%fj = 0, with positive, 75^ = +ip, (right-handed) 
and negative, 75^ = —tp, (left-handed) chiralities respec- 
tively. 

This link between the topological charge and the chi- 
rality (or helicity) of zero eigenmodes of the Dirac op- 
erator identifies parity as the transformation for chang- 
ing the sign of the topological charge while leaving the 
action invariant. Helicity transforms as a pseudoscalar 
under rotations, whereas the action of QCD, designed to 
conserve parity, transforms as a scalar under rotations. 
For the improper rotation of the parity transformation, 
the right-handed modes will become the left-handed, and 
vice versa, changing the sign of Q. 

In deriving the transformation of the links under par- 
ity, we begin with the transformation of the gauge po- 
tential. Under the parity transformation 

Ai(x,t) -» -Ai(-x,t) , and di -> -o\ , (6) 

for the spatial indices i — 1,2,3 and 

A^{x, t) — » A±{—x, t) , and 84 — > 84 . (7) 
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Here the 3x3 colour matrix degrees of freedom of A^x) 
are implicit. The non-abelian field strength tensor F pv is 

F^u = d^A v - d v A ll + i g [A„,A V ] . (8) 

Under parity transformations 



Fij — > Fij , F^j 
Hence the action 
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S=\j d 4 xTv(F^(x)F^(x)) 



(10) 



is invariant under parity while the topological charge 



d 4 x e^p^Tr (F MV (x) F pa (x)) , (11) 
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changes sign due to the presence of e^, vp<r ensuring the 
presence of one and only one time component in the field 
strength tensor product. 

The gauge-field links are related to the gauge potential 
via 

U lt (x)=exp(igJ A^x + X fi) dX^j , (12) 

where x = (x, t) . Under the parity transformation, the 
spatial components, i = 1,2,3, of the links transform as 

VUi(x)T>i = exp(-igj^ A^-x - Xz) dXj ,(13) 
= Uj(-x-ai), (14) 

where 

-x = (-x,t). (15) 
Similarly, the time components transform as 

VU t {x)V ] = expfi$jf A t (-x + Xi)dX^j , (16) 

= U t (-x). (17) 

The exact nature of the parity transformation on the 
links follows from the fact that parity is an exact sym- 
metry on the hypercubic lattice. 

It is interesting to examine the manner in which the 
lattice action density and lattice field-strength tensor 
transform under parity. Consider the product of links 
about an elementary plaquette located at space-time 
point x, in the \i - v plane 

U^(x) ee U^x) U v (x + fx) Ufa + 0) Ut{x) , (18) 

and its transformation under parity. For spatially ori- 
ented plaquettes, the untransformed plaquette Uij(x) 
originates from x, and loops counter-clockwise in the pos- 
itive i, j direction 



U i (x)U j {x + i)U}(x + j)uUx). 



Under the parity transformation Uij(x) transforms to 

Uj(-x-t) u}{-x-i-3) Uii-x-j-i) Uji-x-j) , (20) 

which is the spatial plaquette originating from — x, loop- 
ing in a counter-clockwise orientation again, this time in 
the negative i, j direction. The space-time oriented pla- 
quettes originating from x and looping counter-clockwise 
in the positive i, t direction 



Ui {x) U t {x + i) UJ (x + t) U} (x) 



(21) 



transform to 

U}{-x - i) U t {-x - i) Ui(-x - i + t) U}{-x) , (22) 

which is the space-time oriented plaquette originating 
from —x but this time looping clockwise in the negative 
i, positive t direction. 

The connection between these link products and the 
lattice field strength tensor is provided by 



gF, v {x) = 1 



(23) 



where O flu (x) is the sum of 1 x 1 link paths oriented 
about x in the fi - v plane 

O^(x) = U^x) U v {x + fx) Ufa + 0) Ufa) 
+ U v (x) Ul(x + 0- fi)Ufa- fi) U„(x - fi) 
+ Ufa - fi) Ul(x - fb-0) Ufj,{x - fi - v) U u (x - 0) 
+ Ufa-0)U„(x-0)U v (x + fl-0)Ufa). (24) 

Taking the Hermitian conjugate of O pv (x) changes the 
orientation of the link products from counter-clockwise to 
clockwise. Noting further that F fiu (x) is odd under such 
transformations, we see that the orientations of the link- 
product parity transformations of Eqs. (|20|) and (|22f) are 
precisely those required to recover the continuum trans- 
formations of F l j, l/ (x) in Eq. ©. While the topological 
charge density undergoes a parity transformation, the 
magnitude of the topological charge remains invariant, 
as all lattice sites are summed over. Similarly, the lattice 
action 



5 = / ? E^ Tr ( 2 



U„u(x) - UUx)) 



(25) 



is even under the orientation transformation U pv (x) — > 
UJ lv (x), and since all sites are summed over, the lattice 
gauge action is invariant under the lattice parity trans- 
formation. 



III. 



NUMERICAL RESULTS 



(19) 



The parity transformations of the gauge links in 
Eqs. (jTlJl and lfT7|l have been coded in Fortran 90 [2(j 
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and tested on a 12 3 x 24 lattice with the tadpole-improved 
Liischer-Weisz gauge action [21|. We use three methods 
to examine the properties of the parity transformed con- 
figurations. 

The 3-loop improved field strength tensor [l8| accom- 
panied by the associated 3-loop improved cooling action 
are used to provide a determination of the topological 
charge Q and action S for a cooled gauge configuration 
and its parity transformed partner. The transformed ac- 
tion and topological charge magnitude agree to machine 
precision. 

The low-lying spectral flow of the hcrmitian Wilson- 
Dirac operator H\y — 75-Dw is examined. Under the 
parity transformation, the low-lying eigenvalues of Hw 
change sign such that the slope of the spectral flow (pro- 
portional to the topological charge giving rise to the ze- 
romode) changes sign. Since H\y is of even dimension 
(lattice volume x number of colors x number of Dirac 
indices), det(H w ) = det(7 5 ) x &et(D w ) is unaffected by 
the parity transformation, as required to preserve parity 
in dynamical fermion QCD simulations. 

Finally, the topological charge is determined by count- 
ing the exact zeromodes of the massless overlap Dirac 
operator. Under a parity transformation, the number of 
exact zeromodes with positive (right-handed) and nega- 
tive (left-handed) chiralities exchange as expected. 

IV. DISCUSSION AND CONCLUSIONS 

The manner in which the parity transformation is used 
in practice will depend on the application at hand for 



the gauge fields. Often, exact parity can be enforced 
during the construction of correlation functions by av- 
eraging opposite-sign nontrivial momenta of the lattice 
Green's functions |22|. In other applications, exact par- 
ity can be enforced by doubling the number of gauge field 
configurations, averaging each field configuration with its 
parity transformed partner. Finally, the transform can 
be applied in the Markov chain process itself, applying 
the transform randomly with a 50% probability prior to 
writing the gauge configuration to disk, or equivalently, 
following the reading of the configuration to be evolved, 
from disk. 



Motivated by the lattice index theorem, we have illus- 
trated how a parity transformation on the SU (3) link- 
variables of lattice gauge configurations reverses the sign 
of topological charge while leaving the action invariant. 
The transformation provides an improved unbiased esti- 
mator of the ensemble average and is essential in improv- 
ing the ergodicity of the Markov chain process. 
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